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We perform three-dimensional relativistic hydrodynamical simulations of the coalescence of 
strange stars and explore the possibility to decide on the strange matter hypothesis by means 
of gravitational-wave measurements. Self-binding of strange quark matter and the generally more 
compact stars yield features that clearly distinguish strange star from neutron star mergers, e.g. 
hampering tidal disruption during the plunge of quark stars. Furthermore, instead of forming dilute 
halo structures around the remnant as in the case of neutron star mergers, the coalescence of strange 
stars results in a differentially rotating hypermassive object with a sharp surface layer surrounded 
by a geometrically thin, clumpy high-density strange quark matter disk. We also investigate the 
importance of including nonzero temperature equations of state in neutron star and strange star 
merger simulations. In both cases we find a crucial sensitivity of the dynamics and outcome of the 
coalescence to thermal efi^ects, e.g. the outer remnant structure and the delay time of the dense 
remnant core to black hole collapse depend on the inclusion of nonzero temperature effects. For 
comparing and classifying the gravitational-wave signals, we use a number of characteristic quan- 
tities like the maximum frequency during inspiral or the dominant frequency of oscillations of the 
postmerger remnant. In general, these frequencies are higher for strange star mergers. Only for 
particular choices of the equation of state the frequencies of neutron star and strange star mergers 
are similar. In such cases additional features of the gravitational- wave luminosity spectrum like the 
ratio of energy emitted during the inspiral phase to the energy radiated away in the postmerger stage 
may help to discriminate coalescence events of the different types. If such characteristic quantities 
could be extracted from gravitational-wave signals, for instance with the upcoming gravitational- 
wave detectors, a decision on the strange matter hypothesis and the existence of strange stars should 
be possible. 

PACS numbers; 04.30Tv, 12.39Ba, 21.65Qr 



I. INTRODUCTION 

The equation of state (EoS) of high-density matter and 
therefore the true nature of compact stars has been a 
mystery since the discovery of these objects. Besides the 
possibility of a neutron-proton composition more exotic 
phases have been proposed to appear in the cores of com- 
pact stars (see e.g. [ij for a review, for early considera- 
tions of quark stars see 0,J3|)- The formulation of the 
strange matter hypothesis [j, |5| introduced the possibil- 
ity that compact stars are so-called strange stars (SSs) 
[a, 0| . According to this hypothesis the absolute ground 
state of matter is formed by a quark phase consisting 
of up, down and strange quarks, and compact stars are 
self-bound objects composed of this kind of matter. It 
turns out that SSs are in many ways similar to neutron 
stars (NSs), e.g. in the mass range and compactness, and 
therefore they can be considered as an alternative expla- 
nation for compact stars [l|, . Throughout this article 
the term compact star refers to either a NS or a SS. 

Theoretically as well as observationally, the ques- 
tion whether the strange matter hypothesis is true and 
whether (at least some) NSs are actually SSs, is an open 
issue. Besides the determination of compact star param- 
eters, several possibilities have been proposed to explore 
the observational consequences of the strange matter hy- 
pothesis, including the possibility that small lumps of 
this strange quark matter (SQM) with baryon numbers 



of about lO'^ and more might be abundant in the cosmic 
ray flux (sec e.g. [l3l)- Experiments like AMS-02, which 
is planned to be installed on the International Space Sta- 
tion in 2010, have the potential to detect these so-called 
strangelets [Tl|, [l2j ■ Also, SQM might be detected with 
the LHC at CERN [l3l4Ta|. For a review on these con- 
siderations and other observational implications of the 
strange matter hypothesis see d, [l^ . 

With the current and upcoming GW detectors like 
LIGO [13 and VIRGO [l| and the prospect of a de- 
tection of signals from compact object binaries [l9| , the 
question naturally arises if the signal of mergering NSs 
could be distinguished from SS mergers, especially be- 
cause compact object binaries are considered to be among 
the most promising sources for these detectors. 

The systematic investigation of the imprint of the EoS 
on the gravitational-wave (GW) signal is still in its in- 
fancy. Most studies consider simplified EoSs like poly- 
tropes and try to explore the chances to measure some 
general c omp act star properties like the stellar radius 
(e.g. [20I l2l|). It is not clear whether a decision on the 
strange matter hypothesis could be made on the basis 
of such measurements, because the compactness of SSs 
can be in the same range as that of NSs. First attempts 
of exploring the consequences of microphysical EoSs in- 
cluding nonzero temperature effects in a relativistic treat- 
ment, necessary for reliable results, have been made by 
[2^ . Also fully relativistic studies have been conducted 
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in which temperature corrections were approximated by 
an ideal gas component added to a zero-temperature mi- 
crophysical EoS (e.g. and preceding works). 

Strange stars as sources of GWs have been consid- 
ered in a rotational equilibrium approach investigating 
the (final phase) of the inspiral of SS binaries (24l. [25j. 
Quasiequilibrium orbits were constructed to determine 
the innermost stable circular orbit (ISCO). It was found 
that the frequency of the GWs at the ISCO depends 
on the compactness of the SSs. The orbital evolution 
and the associated GW emission of SS-black hole bina- 
ries has been estimated within a semianalytic model in 
j26| . Moreover, the signals from various instabilities of 
rotating SSs have been worked out in [27l - l29j (see also ref- 
erences therein) . In (30l | it was found that the frequencies 
of the g-mode in newly born SSs are significantly lower 
than those of NSs. Differences of the fundamental pres- 
sure modes are discussed in (sif where a discrimination 
was found only to be possible if additional information 
like the mass was available. Nevertheless it is not clear 
whether these discriminating features are relevant as the 
corresponding GW signals might be too weak for mea- 
surements. 

In this article we examine the GW signals emitted by 
merging SSs. We explore two different EoSs represent- 
ing different possibilities of SQM properties and consider 
the inspiral phase, the final plunge, and the postmerger 
stage. We compare the GW characteristics of SS coa- 
lescence with those of ordinary NSs and discuss if GW 
observations could be decisive on the strange matter hy- 
pothesis. For that we conduct three-dimensional rela- 
tivistic hydrodynamical simulations that allow us to fol- 
low the evolution of the whole merging process and to 
extract observational signatures. To our knowledge this 
is the first study considering these events in dynamical 
simulations; only results for mergers of SS-black hole sys- 
tems have been so far reported in [s^ , where Newtonian 
hydrodynamics were used and the black hole (BH) was 
modeled by a pseudorelativistic potential [ssj . 

The paper is structured as follows: Sect. [ll] will intro- 
duce the underlying model and numerical methods. In 
Sect. lIIII thc properties of the EoSs used in this study will 
be described. The simulations and some general aspects 
will be discussed in Sect.|lVl In Sect. |V] the GW signals 
will be presented and in Sect. IVII we will close with a 
summary and our conclusions. 



II. NUMERICAL METHODS 

For this study we used a version of the general rela- 
tivistic smooth particle hydrodynamics code described in 
[3^ [35} and we refer the reader to these publications for 
details. As in grid-based Eulerian relativistic hydrody- 
namics schemes we introduce a set of so-called conserved 
quantities (p*, Mi, r), the conserved rest mass density, the 
conserved specific momentum and the conserved energy 
density, which relate to the primitive quantities {p,Vi,e), 



the rest mass density, the coordinate velocity and the 
specific internal energy, via 

p* = pau^^j^, (1) 
Ui = hu^ = h{v' + /3')V'^w°, (2) 

Here the Lorentz factor is defined by F = = (1 -I- 
7'''uiMj)"'^/^, and Ui are the time and space components 
of the eigenvelocity, respectively, /i = 1 + Pj p + e rep- 
resents the specific relativistic enthalpy, and the metric 
potentials a, ■0 and /3* are introduced below. 

The evolution of these conserved quantities can be ex- 
pressed in a Lagrangian manner, e.g. comoving with 
the fluid (see [34[[3^). Having at hand this Lagrangian 
formulation of relativistic hydrodynamics the discretiza- 
tion according to the Smoothed Particle Hydrodynamics 
(SPH) method can be done in a straightforward manner 
(see e.g. [13, HI]) by defining a smoothing operator 

{A{r))= J A{r')W{\r-r'\;h)d^r' (4) 

for any quantity A(r). The analytically given normal- 
ized, continuous, differentiablc kernel function W{r; h) 
has compact support and is characterized by the smooth- 
ing length h [33|. Representing the fluid by a sample of 
particles with fixed rest masses the above integral can 
be approximated by a summation 

(A(r))«^yl(r,)^W^(|r-r,l;/i) (5) 

over neighboring particles, where p* is the conserved rest 
mass density of particle "a" located at . Rewriting the 
evolution equations in this form has the advantage that 
partial derivatives can be expressed by means of an in- 
tegration by parts in equation Q as derivatives of W, 
which are analytically given. Therefore one saves the 
computationally expensive evaluations of partial deriva- 
tives. 

In our code the resulting system of ordinary differen- 
tial equations (see HH) is solved with a fourth order 
Runge-Kutta method. In every timestep one recovers 
the primitive quantities (p, Vi, e) from the evolved quanti- 
ties {p*,Ui,T) by solving numerically the system ((T])-([3]), 
where the temperature T is determined from the spe- 
cific internal energy e by means of the EoS. The EoS 
P{p,Ye,T) finally closes the system of hydrodynamical 
equations and will be discussed in Sect. IIIIl The initial 
electron fraction of the cold stars is defined by the con- 
dition of beta-equilibrium and is assumed to be advected 
with the fluid (^^ = 0) because the timescale of lepton 
number transport (at least for the bulk of the stellar mat- 
ter in the interior, where the GW signal is determined) 
is long compared to the dynamical timescale. This, how- 
ever, does not hold for the matter shed off the surface of 
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the merging stars, where lepton number and energy loss 
may be important. 

In contrast to standard Newtonian SPH, gravity can- 
not be implemented within the particle picture (e.g. with 
a tree method). Instead, the Einstein field equations have 
to be solved. For this we employ the conformal flatness 
approximation, which assumes that the metric can be 
written as ds^ = (— + I3if3^)dt^ + 2(3idx^dt + y^^dx^dx^ 
with 7ij = -0^(5^. Here a is the lapse function, /?' is the 
shift vector and 7ij is the spatial part of the metric ap- 
proximated by the conformal factor ip and the Kronecker 
delta ^ij [H, ll^ . Within this approach the Einstein field 
equations reduce to a set of five coupled nonlinear ellip- 
tic partial differential equations with noncompact source 
terms. These equations are discretized on a grid covering 
a domain around the binary, respectively the postmerger 
remnant, and are solved with a full multigrid method. 

Since GWs are neglected by this ansatz for the metric, 
a scheme simulating the backreaction of the GW emis- 
sion has to be implemented in order to account for the 
loss of energy and angular momentum carried away by 
the gravitational radiation, which drives the binar y to in- 
spiral and to finally merge. Following the ideas of pol.l4]| 
this is done in a post-Newtonian framework. For details 
of the implementation of the GW backreaction see [s^ , 
where the code we use is described in all parts, including 
the treatment of shocks. 

Magnetic fields are not included in our implementa- 
tion. The effects of magnetic fields smaller than lO^^G 
on the merger dynamics are negligible on shorter than 
secular timescales 142, H . 



A. Initial data 

Prior to the actual merger simulations the binaries are 
evolved without the effects of the GW backreaction in 
order to find the corresponding quasiequilibrium state of 
the configuration. This is achieved by the relaxation of 
the particle distribution, which is slightly off equilibrium 
after setting up the two individually hydrostatic stars in 
a close orbit at the beginning. 

The initial velocity profile is assumed to be irrotational 
because the shear viscosity of SQM is comparable to that 



of nuclear matter [4J] . Therefore the time during inspiral 
is not sufficient to yield tidally locked systems [43. Also 
because of time arguments it is reasonable to start with 
initially cold SSs in beta-equilibrium. 



III. EQUATION OF STATE 

The EoSs of SQM are derived within the MIT bag 
model 0, Initially developed to describe nucleons, 
the application to bulk quark matter is rather simple. In 
a certain volume, the bag, the quarks are considered to be 
free Fermi gases. The masses of the up and down quarks 
can be neglected, while for the strange quark mass a value 



of ras — 100 MeV is adopted |48|. The confinement of 
QCD is simulated by a nonzero pressure of the vacuum 
outside onto the bag given by the bag constant B. Within 
the MIT bag model the pressure P, the energy density e 
and the baryon density p read @ 



+7l(fc, -/i/)] - S, 



(6) 



E 



/ 



If 

2^^ 7o 



dk k'^ e -f {k)[n{k , iij) 



+nik,^tif)]+B, (7) 

where e/(fc) = {mj + fc^)^/^ and the Fermi-Dirac distri- 
bution function n{k, Pf) = (exp ([e/(fc) — Pf] /T) + 1)^^ 
for a temperature T is used. The chemical potentials of 
the different quark flavors are denoted hy fif. Note that 
we set c = h = 1. The sum extends over all quark flavors 
with the degree of freedom 7/ = 2 x 3 for two spin states 
and three color states. These equations can be solved 
analytically in two limits: for massive quarks and zero 
temperature, or alternatively for massless quarks with 
nonzero temperature Q. However, in this study a mas- 
sive strange quark is assumed and thermal effects are sup- 
posed to play a role during the merging process of SSs and 
therefore the above set of equations needs to be solved 
numerically j49j . The conservation of baryon charge and 
electric charge and the condition of electric charge neu- 
trality lead to a set of constraints for the chemical poten- 
tials under which the above equations are solved. These 
constraints also determine the chemical potential of the 
electrons, which are treated independently as an ideal 
Fermi gas. For details we refer to jsl. Isoj. 

The bag constant is a free but constrainted parameter 
of the model. Because a spontaneous conversion of nu- 
clei into SQM is not observed, the lower limit is given 
by i? = 57 MeV/fm"^. In order to obtain an energy per 
baryon lower than the one of ordinary nuclear matter the 
upper bound of B equals 84 MeV/fm'^. The range of B 
can slightly differ for other choices of ms or if one in- 
cludes interactions between the quarks. These limits for 
B correspond to an energy per baryon lower than that 
of nuclear matter {E/A = 930 MeV) and therefore they 
account for absolutely stable quark matter, i.e. the true 
ground state of matter. In this study we use the con- 
stant values B = 60 MeV/fm^ (E/A ^ 860 MeV) and 
B = 80 MeV/fm3 {E/A = 921 MeV) and refer to these 
EoSs as MIT60 and MIT80. Because of these choices 
of B the EoSs roughly represent the two extreme cases 
of the underlying microphysical model. Potential exten- 
sions of the model could for instance account for the ef- 
fects of color superconductivity or interactions between 
the quarks. In order to have a clear parameter depen- 
dence we omitted the inclusion of such corrections, which 
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anyway could be absorbed in an "cfFective bag constant" 

In principle SSs can carry a thin nuclear crust whose 
maximum density cannot exceed the neutron drip density 
Q- When free neutrons form they settle gravitationally 
and are converted into SQM when they come in contact 
with the quark phase. For this reason the mass of the 
crust is only of the order of 10~^ Mq. Therefore we ne- 
glect it in our simulations, because it is not of importance 
for the dynamics of the system. 

Since we want to compare our findings for SS mergers 
with NS coalescence, we also performed simulations with 
EoSs describing NS matter. To this end we employed the 
same nuclear EoS as used in [s^, namely those of Refs. 
[5^ [5^ (see there for more details) and abbreviate them 
in the following by Shen and LS, respectively. These 
models also span a considerable range of variation for NS 
matter in the sense that LS yields fairly compact stars 
while Shen results in significantly less compact objects 
(Fig. El). 

Figure [1] shows for all EoSs used in this study the pres- 
sure as a function of the rest mass density p = rriunB {ns 
is the baryon number density, mu the atomic unit mass) 
for representative temperatures of 0.1 MeV, 10 MeV and 
30 MeV. For the EoSs of SQM one can clearly see the 
drop of the pressure towards zero at densities well above 
nuclear density (2.7 • 10^** g/cvn?). Moreover, in the case 
of SQM, except for a narrow range of densities close to 
this minimum density, the temperature has only a minor 
infiuence on the pressure. This is due to the fact that 
SQM is forming a degenerate Fermi gas of ultrarelativis- 
tic particles. Also for the nuclear EoSs thermal effects 
become more important the lower the density is. At a 
density of 10^^ g/cm^ the increase of the temperature 
from 0.1 MeV to 30 MeV leads to a pressure increase of 
about one order of magnitude. 

Solving the Tolman-Oppenheimer-VolkofF stellar struc- 
ture equations for the EoSs introduced above, i.e. solving 
the general relativistic equations of hydrostatic equilib- 
rium, one obtains the mass-radius relations of bare SSs 
in beta-equilibrium at T = 0. Figure [5] shows the cor- 
responding results with M referring to the gravitational 
mass (Arnowitt-Deser-Misner mass) in isolation and R 
being the radius in Schwarzschild coordinates. First of 
all, one recognizes the typical inverse mass-radius rela- 
tion for the NS models, which is much different from 
that of SSs. The maximum mass of MIT60 and MIT80 
is 1.88 Mq and 1.64 Mq, respectively. Remarkably, the 
maximum mass of LS is 1.83 A/q, thus comparable to the 
one of MIT60, while the Shen EoS supports stars with 
masses up to 2.24 Mq. More massive nonrotating ob- 
jects collapse to BHs. Note that the inclusion of rotation 
increases these maximum masses. Rigidly rotating stars 
more massive than these maximum nonrotating objects 
are called supermassive, while differential rotation, which 
increases the maximum mass further in comparison to 
uniform rotation, allows for so-called hypermassive stars 
exceeding the maximum mass of supermassive stars [s^ ] . 
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FIG. 1: Pressure as a function of rest-mass density p = muTiB 
for the MIT60, MIT80, LS and Shen EoSs, and for different 
chosen temperatures. The EoSs for SQM can be easily recog- 
nized by tlieir rapid decrease of the pressure below a certain 
density. The electron fraction for the nuclear EoSs is set to 
0.05 and for the quark matter EoSs to 2.5 ■ 10~^. 



The stellar radii for LS and MIT60 in the high-mass range 
are similar and therefore the compactness GM/c^R, too. 
In general, SSs are more compact than NSs, and MIT80 
yields smaller stars than MIT60. 

In order to illustrate the effects of nonzero tempera- 
ture, we also plot the mass-radius relations for isentropic 
stars with nonzero entropy. The entropy is fixed by the 
value obtained when matter at about three times nuclear 
saturation density is chosen to have a temperature of 30 
MeV, which are typical temperature and density values 
reached during a binary merger. This corresponds to an 
entropy per baryon of 1.1 for the Shen EoS and 1.5 ks 
for the LS EoS. The same criterion gives an entropy per 
baryon of 3.2 ks for the MIT60 EoS and 3.1 ks with the 
MIT80 EoS. One can see that thermal effects increase 
the maximum mass of NSs only slightly. However, the 
influence on the NS radius is bigger the smaller the mass 
is. Especially at lower densities thermal effects become 
important as mentioned above. Remarkably, doing the 
same for the SS EoSs yields mass-radius curves that are 
hardly distinguishable from those for T = 0, which is also 
understandable from Fig. [1] 

We stress the fact that SQM as a small nugget but 
also as an SS is self-bound by QCD interactions (and a 
disintegration of the constituents is energetically not fa- 
vorable). As the pressure goes to zero the density at the 
surface of an SS (or a small nugget of SQM) adopts a 
value of the order of the nuclear saturation density (see 
Fig.[T]). It implies that this type of matter cannot become 
arbitrarily dilute. This also means that for the binding of 
a SS gravity is not essential, it only compresses the star 
additionally. In contrast, an ordinary NS is solely bound 
by gravitation and would "explode" without the effect of 
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FIG. 2: The stable branches of the mass-radius relations for 
compact stars with EoSs MIT60, MIT80, LS and Shen. M 
refers to the gravitational mass in isolation and R is the ra- 
dius in Schwarzschild coordinates. The black curves display 
the mass-radius relations for T = Q, whereas the red curves 
show them for isentropic stars with an entropy per nucleon 
of 1.1 fcs for the Shen EoS (right curve) and 1.5 fcs for the 
LS EoS (left curve). For the SQM EoSs the differences due 
to thermal effects are significantly less than 1% even for the 
chosen entropy of 3.2 fcs in the cases of the MIT60 EoS and 
3.1 fcs for the MIT80 EoS (see text for an explanation of the 
choices of the entropy values) . 



gravity. Pictorially spoken, SQM behaves hydrodynami- 
cally like a liquid with surface tension. 

Finally, we note that the mass-radius relation is a cru- 
cial quantity determining the dynamics and the outcome 
of a merger (see below and [sBI, fssl)- In our SQM model 
the bag constant is the only parameter varied and thus af- 
fects entirely the mass-radius relation. Therefore our re- 
sults can be interpolated for other choices of the bag con- 
stant between the two considered nearly extreme cases. 
Similarly, one can estimate how disregarded effects like 
quark interactions and color superconductivity might in- 
fluence our results by inspecting their effects on the mass- 
radius relation (see also [1^). In the same way one can 
judge on the impact of other quark matter models be- 
yond the standard MIT bag model. In Ref. [13, for 
instance, running quark masses were considered, which 
affect the mass-radius relation such that SSs can have 
a higher maximum mass, which will lead to differences 
compared to our models with respect to the onset of BH 
formation. 



IV. SIMULATIONS 

The mergers of SSs fall into two categories similar to 
the case of NS coalescence. For high total binary masses 
the remnant of the merging stars cannot be supported 
against gravitational collapse and a BH forms shortly af- 
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FIG. 3: Time evolution of the maximal rest mass density (up- 
per panel) and the minimal lapse function (lower panel). Note 
that we added a time shift of -4.5 ms to the 1.2Af0-|-1.35M0 
model and of -6 ms to the others to avoid overlapping curves. 



ter the stars have come into contact. The forthcoming 
appearance of a BH is indicated by a steep increase of 
the central density pmax exceeding more than twice the 
maximum central density of nonrotating SSs accompa- 
nied by a steep decrease of the lapse function a, both on 
a timescale smaller than the sonic time scale. This can 
be seen in Fig. [3] (dotted line) for a symmetric binary 
with two 1.8 Mq stars described by the MIT60 EoS. On 
the other hand, a merged object with a sufficiently low 
total mass can be temporarily stabilized mainly by dif- 
ferential rotation; a corresponding hypermassive object 
forms when the total mass exceeds the maximum mass 
of supermassive stars [l^l- This object will also collapse 
to a BH after angular momentum redistribution. 

Figure |3] shows the different stages of an SS merger, 
which finally leads to the formation of such a hypermas- 
sive object. The stars have a mass of 1.2 Mq and 1.35 Mq 
and they are described by the MIT60 EoS. The orbital 
period of the counterclockwise rotating binary shortly 
before the merger is of the order of 2 ms (upper left 
panel). Subsequent loss of angular momentum and en- 
ergy due to GW emission during the inspiral phase leads 
to a shrinking orbit until the stars finally merge with 
a relatively big impact parameter (upper right panel). 
An initially deformed but finally approximately axisym- 
mctric differentially rotating hypermassive object forms, 
which is stable for more than 10 ms (lower panels) (for a 
discussion of the delay time Tdciay to BH collapse see [ssj 
and Sec. IVD[) . After several revolutions of the merger 
remnant tidal arms form and matter from the tips of 
these tails becomes gravitationally unbound and so con- 
tributes to the cosmic ray flux of strangelets. For the 
implications of this aspect we refer to [5^ • Most of the 
matter in these spiral arms ends up in orbits around the 
remnant and forms a geometrically thin accretion disk 
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(see Figs, m and [6]). 

One can also sec that the matter that is not part of 
the hypermassive object forms clumps, which is a conse- 
quence of the fact that SQM is self-bound. NS mergers 
instead form dense remnants surrounded by inflated halos 
with toroidal structures and a smooth density distribu- 
tion (see Figs. Eland El). 

Since we want to discuss in this paper how to discrimi- 
nate between SS mergers and NS mergers we will concen- 
trate in the following on a selection of models for the LS 
and MIT60 EoSs and consider general features to reveal 
fundamental differences. The reason to focus on these 
two EoSs is that LS as nuclear EoS and MIT60 as SQM 
EoS are very similar at densities around 10^^ g/cw? and 
lead to compact stars with similar radii for Af > 1.5 A/© 
(see Figs. [Hand El. The Shen EoS and the MIT80 EoS 
exhibit bigger differences and therefore one expects that 
the distinguishing features are more pronounced. In to- 
tal, we have computed more then 50 models and discuss 
in more detail 6 simulations of the 9 listed with their 
properties in Table HI The choice of three of these binary 
configurations is motivated by the fact that population 
synthesis studies (soj predict a "mean" binary with a to- 
tal mass of about 2.7 Mq and a mass ratio close to unity, 
which is consistent with the masses of the well measured 
compact star binaries (60j . In addition, we include asym- 
metric systems with 1.2 Af© and 1.35 Mq to investigate 
the effects of mass ratios unequal to unity. We choose 
relatively low masses because only in these cases hyper- 
massive objects can form. Finally, we consider a binary 
configuration in which the two components have approx- 
imately the same radius for the LS and MIT60 EoSs, i.e. 
we take the stellar masses to be close to the maximum 
mass for the LS and MIT60 EoSs (see Fig. We wih 
discuss for these models the general dynamics and out- 
come of the merging phase, the ejecta and torus masses, 
the influence of thermal effects, and in Sect. [Vlthe fea- 
tures of the GW signals. In particular in the cases where 
a hypermassive object forms, the GW emission from the 
postmerger phase yields characteristic information about 
the EoS of matter at very high densities. 

In Table U we also list results for the Shen EoS without 
discussing these models in much detail, because the merg- 
ers described by this EoS behave qualitatively similar to 
the ones with the LS EoS. The comparison between dif- 
ferent EoSs of NS matter was reported in [sHi • We refrain 
from presenting results for the MIT80 EoS, because all 
binary configurations considered in Table U form a BH 
promptly when computed with the MIT80 EoS. How- 
ever, we will address results with this EoS extensively in 
Sect. El 



A. Binaries with 1.2 Mq and 1.35 Mq 

The two 1.2A'/Q-t-1.35A'/Q binaries with the LS and 
MIT60 EoS both form hypermassive objects. The general 
dynamics of the SS model were described above. FigureEl 



shows several stages during the merging of the NS binary 
to be compared with Fig SI Different from the SS case, 
the less massive star is stretched during the final stage of 
the coalescence and directly forms a massive spiral arm 
(upper right panel), which means that the lighter star is 
tidally disrupted during the merging. Continuous mass 
shedding from this spiral arm feeds an extended torus- 
like halo around the remnant (lower left panel). The 
quark stars in contrast, coalesce as entire stars without 
the prompt formation of such a feature (see Fig. [71). 

The structure of the surface layers of the merging stars 
and of their remnants is fundamentally different. In the 
case of the LS EoS, matter is also squeezed out from the 
contact layer, which does not happen for MIT60. Ad- 
ditionally, shortly after the first contact matter streams 
off from the whole surface of the remnant and fills the 
diffuse halo (Fig. [HI). 

In contrast to the NS case, during the complete merg- 
ing process of the SSs first the stars then the merged 
object retain their sharp boundary, no matter is spread 
into the surroundings. Recall the liquidlike behavior of 
SQM. Only relatively late after about two revolutions gas 
is stripped off from the remnant due to the formation of 
two filigree but very dense spiral arms. This process is 
constrained to the close vicinity of the equatorial plane 
(see Fig.[6l) and subsequent fragmentation leads to a disk 
consisting of clumps of SQM orbiting the compact rem- 
nant. Only a small amount of matter from the tips of 
the tidal tails can escape from the merger site. Remark- 
ably, we find a gap between the merger remnant and the 
disklike structure (which due to the projection of the par- 
ticle positions onto the x-z plane is not visible in Fig.[ni). 
To our knowledge such a distinct feature is unique to SS 
mergers. The differences of the coalescence and the evolu- 
tion of the merger remnant can be understood by means 
of the higher compactness and self-binding of SQM. We 
suspect that the reason for the gap is the occurrence of 
an ISCO outside the surface of the central remnant. An 
IS CO larger than the stellar radius is known as generic for 
slowly as well as rapidly uniformly rotatingSSs, in partic- 
ular for supermassive configurations |63| . Therefore, 
we expect that this holds also for hypermassive differ- 
entially rotating objects. Note that the omission of a 
possible nuclear crust of the SSs does not change this 
picture because the crust would contain only very little 
mass. A halo fed by this material would be much more 
dilute than the one in the NS case. 

In addition, with the LS EoS the central lapse func- 
tion decreases during the ringdown of the hypermassive 
object, while with MIT60 it reaches quickly a nearly con- 
stant value as can be seen in Fig. [31 This behavior is 
compatible with the central density, which increases for 
the remnant consisting of nuclear matter described by 
the LS EoS (see Fig. [31) . 

Two quantities can be extracted from our simulations 
for characterizing the mass shedding during the merging: 
One is the amount of matter that becomes gravitation- 
ally unbound, the other is the mass of a torus that will 
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TABLE I: Models that are discussed in the text. Mi and M2 refer to the individual gravitational masses in isolation of the 
two components of the binary, /max, /peak and /gap are characteristic frequencies of the GW signal (see text for definitions). 
Tdeiay is the delay time between the merging and the formation of a BH. Note that in the cases where the BH collapse does 
not occur during the simulation, we give a lower limit for Tdday, this value is determined by the finite simulation time and 
may differ significantly from the true value (see [HH]). Mtorus and Mejccta are the torus and ejecta masses at the end of the 
simulations. Note that the torus masses of models LS I.SSMq + I.SSM© and Shen 1.35M0 + 1.35M0 are still increasing at the 
end of the simulations. In the cases we find a vanishing torus mass or ejecta mass we provide the upper bound of this quantity 
given by the mass resolution of our simulations. Tmax denotes the maximal temperature reached during the evolution. For the 
I.8M0 + I.8M0 configurations the presented temperatures are given by the last dump step during the simulations. The actual 
maximal temperatures during the merging process shortly before the horizon formation are expected to be higher. Table entries 
in parentheses correspond to models where zero temperature was imposed. 
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remain after the remnant has collapsed to a BH. The de- 
termination of these two quantities is described in (ssj . 
They have direct physical relevance because the ejecta 
mass is important for nucleosynthesis contributions from 
NS mergers and for strangelct injection to cosmic rays in 
the case of SS mergers as discussed in [HBl • On the other 
hand, the torus is a potential energy reservoir for power- 
ing gamma-ray bursts. The implications of a SQM disk 
or torus for explaining gamma-ray bursts has not been 
explored yet. A discussion of this question is beyond the 
scope of the present paper because it requires an under- 
standing of the neutrino cooling and viscous evolution of 
the SQM disk. 

We find a torus mass of 0.14 Mq for LS and 0.08 Mq 
for MIT60 (see Table IJ). The differences can be under- 
stood by the different dynamical behavior of the merg- 
ing stars, where tidal disruption favors the formation of 
a more massive torus. For a detailed discussion of the 
torus masses and their dependence on the NS merger dy- 
namics we refer the reader to To clarify this point 
we show the evolution of the torus mass for both EoSs 
in Fig. [T] As one can see by the arrows indicating the 
merger time, the rapid increase for the LS EoS is due 
to the tidal disruption of the lighter star (compare with 
Fig. [S]). In contrast, the torus mass of the quark star 
remnant increases relatively late in the evolution by re- 
distribution of angular momentum to the outer parts of 
the remnant (compare with Fig. 2]). 

The ejecta masses are 0.008 Mq for LS and 0.004 Mq 
for MIT60 (Table m. The smaller values in the latter 
case can be explained as a consequence of the different 
dynamical behavior due to the higher compactness and 



the self-binding of the SSs, also impeding the ejection of 
SQM from the contact layer. This is confirmed by Fig. [5] 
showing the time evolution of the ejecta masses. While 
during the NS merging matter becomes unbound shortly 
after the first contact (compare with Fig. [5]), in the SS 
case the ejecta mass starts growing only about 2 ms after 
the plunge when the spiral arms form (Figs. U] and [T]). 
Note that not all SPH particles that formally fulfill the 
ejecta criterion (dashed black curve in Fig. |5]) are finally 
able to escape to infinity because they interact with other 
particles. Here again the "stickiness" of SQM as self- 
bound fluid manifests itself. By applying an additional 
distance criterion we estimate the true amount of ejecta 
(solid black curve in Fig. [5]) 



B. Binaries with 1.35 M© and 1.35 Mq 

The merging of two stars with 1.35 M© and 1.35 Mq 
proceeds differently from the asymmetric case for the nu- 
clear LS EoS. In contrast, in the case of quark-star bina- 
ries the merging of symmetric and asymmetric systems 
is similar, i.e. the stars collide as whole objects with 
only slight deformation (no tidal disruption). Since the 
stars possess the same mass, the coalescence and rem- 
nant dynamics are symmetric with respect to the two 
stars. A double-core structure forms and the cores os- 
cillate against each other until the remnant settles to an 
axisymmetric hypermassive object. However, after about 
5 revolutions and several bounces of the two high-density 
cores, the SS merger remnant collapses to a BH. The time 
Tdciay between the plunge (time when the GW amplitude 
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becomes maximal) and the collapse is 4.37 ms. The for- 
mation of spiral arms begins after about 2.5 revolutions 
and ends shortly before the collapse. 

Again the way matter is shed off from the merged ob- 
ject is remarkably different. While in the SS case promi- 
nent, fine spiral arms develop that are constrained to the 
vicinity of the equatorial plane, in the NS merger case 
matter is stripped from the whole surface of the rem- 
nant. The spiral arms of the SS merger remnant appear 
again relatively late and only via these tidal tails matter 
is carried away from the merger remnant before the lat- 
ter collapses to a BH. In contrast, the merging NS binary 
does not develop pronounced spiral arms. 

The estimated torus and ejecta masses are given in 
Table H The torus mass with the LS EoS is significantly 
lower (0.04 Mq) than the torus mass with the MIT60 
EoS, but it is still increasing at the end of the simulation. 
The torus mass of the SS merger is not reliable because 
the collapse to a BH occurs in a phase where stationarity 
has not been achieved. The value should be considered 



only as a rough estimate. 

The amount of matter that becomes gravitationally 
unbound from the merger site is 0.002 Mq for the NS co- 
alescence and 0.001 Mq in the case of the SS binary. The 
differences are plausible in view of to the higher compact- 
ness of SSs and the special properties of SQM, where ad- 
ditional energy is required to overcome the self-binding. 
In comparison to the asymmetric binary configurations 
discussed above, the torus masses as well as the ejecta 
masses are lower for symmetric mergers. 



C. Binaries with 1.8 Mq and 1.8 A4q 

The last binary configuration that we include in this 
detailed discussion consists of two stars with the same 
gravitational mass of 1.8 Mq and the same radius for 
both EoSs (about 10.5 km). Both merging binaries 
promptly collapse into a BH after the stars have come 
into contact (sec Fig. [3] for the model with the MIT60 
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EoS). The delay times for that to happen are of the or- 
der of half a millisecond (Table |l| . The stars approach 
each other symmetrically and they become only slightly 
deformed before they touch. Note the central decompres- 
sion phase as reaction to the tidal fields during the merg- 
ing shortly before the central density increases rapidly 
(shown for the MIT60 EoS in Fig. Also if we can- 
not follow the complete collapse to a BH, we are able to 
give an estimate of the torus masses (see [s^ ) . For both 
EoSs we find a vanishing torus mass. This agreement is 
reasonable because the mass and radius of the compact 
objects with both EoSs are the same and the stars are 
only slightly deformed prior to the collapse. With the LS 
EoS we find about 0.0001 Mq that fulfill the criterion for 
gravitationally unbound matter, but we stress that this 
number is very uncertain. These ejecta originate from 
the contact layer between the two stars. In contrast, we 
do not find any ejecta for MIT60 as during a SS coales- 
cence no matter is squeezed out from the region between 



the stars. Strange matter becomes only unbound from 
the tips of tidal tails that form only after several revo- 
lutions of the merger remnant. In the case of a prompt 
collapse, however, this origin is excluded. 



D. Thermal effects 

It is known that thermal effects play a role during the 
merging and postmerging evolution of binary NSs [ssj . 
To which extent this is true in the case of SSs is a priori 
not clear because of the different properties of SQM in 
comparison to NS matter (see Fig. [Ij. Temperature ef- 
fects clearly become important in the contact layer of 
the colliding stars, where compression is extreme and 
shocks occur. Because of the differential rotation the 
hot sheared matter is spread over the entire star forming 
an inhomogeneous temperature distribution in the rem- 
nant. For the MIT60 EoS the maximum temperature 
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FIG. 6: SPH-particle distributions projected onto the x-z-plane perpendicular to the equatorial plane. The left panel shows 
the SS merger remnant of a 1.2M0 + 1.35Mq binary described by the MIT60 EoS, the right panel the corresponding NS merger 
remnant described by the LS EoS. Rest mass densities are ascribed to the individual particles and color-coded, but note that 
due to projection effects foreground particles dominate the visible density structure in the core of the merger remnants. Note 
the logarithmic scale of the rest mass density for the NS merger remnant. The displayed times correspond to the last panels in 
Fig. |4] and Fig. (5] respectively. The plots were created with the visualization tool SPLASH [Glj . 
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FIG. 7: Time evolution of the amount of matter fulfilling the 
torus criterion for the 1.2M0-I-1.35M0 binary mergers with 
the MIT60 EoS (black) and LS EoS (red). The time labels 
of Figs. |4] and [5] correspond to the time measurement on the 
horizontal axis. The arrows mark the merging time defined 
as the moment when the GW amplitude becomes maximal 
(black for the MIT60 EoS; red for the LS EoS). 



reached during the merging of the binary with 1.2 Mq 
and 1.35 Mq stars is about 57 MeV, while the symmetric 
binary with two stars of 1.35 reaches a maximum tem- 
perature of about 65 MeV. In comparison, NS mergers 
described by the LS EoS yield maximum temperatures of 
about 93 MeV and 104 MeV for the same binary config- 
urations. For both kinds of matter the highest temper- 
atures are obtained shortly after the first contact of the 



FIG. 8: Time evolution of the ejecta masses for the 
1.2M0-I-1.35M0 binary mergers with the MIT60 EoS (black) 
and LS EoS (red). The time labels of Figs.|4]and[5]correspond 
to the time measurement on the horizontal axis. In the case 
of the SS merger some SPH particles fulfill the ejecta criterion 
without being able to leave the merger site. The dashed line 
gives the total amount of matter fulfilling formally the ejecta 
condition. The amount of SQM that finally escapes from the 
merger site is shown by the solid line, where we apply as ad- 
ditional criterion that the SPH particles have a distance of 
more than 370 km from the remnant. The arrows mark the 
merging time defined as the moment when the GW amplitude 
becomes maximal (black for the MIT60 EoS, red for the LS 
EoS). 
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FIG. 9: Rest mass density perpendicular to the orbital plane of a merging NS binary with 1.2 Mq and 1.35 Mq components 
for the LS EoS. The left panel shows the results of a simulation where a zero temperature constraint was imposed. The right 
panel gives the rest mass density of a simulation including thermal effects and corresponds to the lower right panel in Fig. [5] 
and the right panel in Fig. [G] The plots were created with the visualization tool SPLASH [gj] . 



stars, when the maximal densities peak for the first time 
(see Fig. [3]). The higher values in the models with the 
LS EoS can be understood by the shock heating and the 
compression of low density material from the surfaces of 
the initial stars. The surface layers of SSs have already 
initially supernuclear densities and therefore their rela- 
tive compression is much smaller. 

In order to check the influence of nonzero tempera- 
tures we performed simulations setting T = for the 
1.2Mq+1.35Mq binaries and the 1.35MQ-hl.35Mo bi- 
naries. This choice corresponds to the extreme and un- 
realistic case of a perfectly efficient cooling that carries 
away all thermal energy instantaneously. A similar ap- 
proach was chosen in [65| with a poly tropic EoS, because 
effectively, the use of P = Kp^ with K,T = const im- 
plies the same physical assumption ("isentropic case"), 
which was compared to results obtained with an EoS of 
the type P = pe including thermal effects ( "ideal- fiuid 
case" ) . In ReL [6^ it was reported that due to the larger 
pressure support in the ideal-fluid case nonzero temper- 
ature leads to a later BH formation than the isentropic 
treatment. Moreover they found remarkable differences 
of the remaining tori after the BH collapse had occurred. 
In particular, in the cases of neglected thermal effects the 
tori were flat, disklike structures, whereas the mergers de- 
scribed by an ideal-fluid EoS produced inflated tori. For 
NS models these flndings are conflrmcd by our simula- 
tions as displayed in Fig. |9] showing the rest mass density 
in the x-z-plane perpendicular to the orbital plane for the 
1.2M0-I-1.35M0 merger remnant with the LS EoS. In the 
left panel T = was imposed, leading to a much smaller 
vertical extension of the outer remnant parts. Note also 
that the density in the disk in the equatorial plane is 
higher by one order of magnitude compared to the den- 
sity in the "hot" torus (right panel). This leads to the 



conclusion that the additional pressure support by ther- 
mal effects is essential for the structure of the remnant. 
We find about 30% higher torus masses in our T = 
simulation with the LS EoS. We attribute this finding at 
least partially to the decrease of the gravitational mass 
due to "cooling losses", which is about 0.035 Mq for 
this configuration. Lowering artificially the gravitational 
mass in our torus criterion by this amount of matter in- 
creases the estimate for the torus mass in our model with 
thermal effects. Qualitatively the same behavior is found 
for the 1.35M0-f 1.35AfQ binary with an increase of the 
torus mass of about 30% in the T = simulation com- 
pared to the T ^ case. 

Another interesting effect of the perfectly efficient cool- 
ing can be seen in the simulation of the binary with two 
1.35 Mq SSs for the MIT60 EoS. As described above, the 
merger remnant of this configuration collapses to a BH 
after Tdciay ~4 ms. Choosing the same binary setup and 
the same EoS but imposing the zero temperature con- 
straint during the evolution, the collapse does not occur 
within the simulation. The high sensitivity of Tdoiay on 
the total mass of the system (see also [11] ) suggests that 
this finding can be explained by the lower gravitational 
mass of the "cold" model. The difference in the gravita- 
tional mass is about 0.06 Mq when the remnant of the 
"warm" simulation collapses. Thus, in contrast to NS 
mergers, where a longer delay timescale for the BH for- 
mation was obtained when thermal effects were included 
j65| . we see the opposite effect in simulations with SQM. 
This can be understood from the fact that thermal energy 
adds significantly to the gravitational mass, but hardly 
affects the pressure and the structure of SSs and their 
merger remnant (see Fig. [5]), because the objects are so 
dense that typical temperatures do not have a big impact 
(see Fig.H]). 
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These conclusions are confirmed by the structure of the 
hot merger remnants described by the MIT60 EoS. The 
hypcrmassive object of the 1.2M0+1.35Mq model (see 
Fig. m and left panel of Fig. ^ looks very similar to the 
cold one (not shown) and qualitatively similar to the re- 
sulting structure of the zero temperature remnants with 
the nuclear LS EoS (left panel of Fig. |9|). The central 
cores of the cold and the warm SS merger remnant both 
have an equatorial diameter of about 20 km and a verti- 
cal diameter of about 10 km, which is in accordance with 
Fig. [21 where warm SQM hardly changes the mass-radius 
relation of SSs. Both mergers end up with a vertically 
fiat disk structure similar to what is found for models of 
NS coalescence with thermal effects neglected (left panel 
of Fig. [9]). As in the NS merger simulations the estimated 
torus masses arc higher for the zero temperature simula- 
tions of merging SS (about 100% for the 1.2A'f04-1.35M0 
model; the comparison for the 1.35M0-l-1.35Af0 model 
is not meaningful because of the uncertainties associated 
with the early collapse of the warm model). 



V. GRAVITATIONAL WAVES 
A. Characteristic features 

In Fig. [inland Fig. [H the GW amplitudes of the plus 
polarization for SS mergers (upper panels) and for NS 
mergers (lower panels) are shown as measured perpen- 
dicular to the orbital plane at a distance D of 20 Mpc. 
The signals were computed for the models discussed in 
detail above, thus for mergers of binaries with 1.2 Mq 
and 1.35 Mq stars and for mergers of 1.35M0-I-1.35M0 
binaries, both with the MIT60 EoS and the LS EoS. They 
can be considered as typical of the cases when a hyper- 
massive object forms. When collapse to a BH sets in, 
our simulation needs to be stopped and the GW signal 
cannot be tracked to its completion. However, simula- 
tions of NS mergers that can follow the formation of the 
BH reveal that the wave signal is strongly damped in 
comparison to the hypcrmassive object case. The only 
remaining signature after the collapse is associated with 
the ringdown of the BH, which is fundamentally differ- 
ent, i.e. it has higher frequencies, lower amplitudes and 
shorter damping times [2J, [gHI • 

The waveforms in Figs. HOl and [121 are extracted by 
an expression of the quadrupole formula that takes into 
account post-Newtonian effects [40j . 

The inspiral phase showing an increase of the frequency 
and of the amplitude can be clearly identified. Its basic 
signal properties can be easily understood by an analyt- 
ical model of orbiting point masses. The corresponding 
wave train of the cross polarization looks similar but is 
phase-shifted by 7r/2. A characteristic quantity that can 
be extracted from the inspiral phase is the maximum 
frequency at the end of this stage just before the merg- 
ing of the two compact stars. More precisely, we define 
/max as the frequency at the moment when the amplitude 
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FIG. 10: Gravitational-wave amplitudes for the plus polar- 
ization measured perpendicular to the orbital plane at a dis- 
tance of 20 Mpc for binary mergers with Mi — 1.2 Mq and 
M2 = 1.35 Mq, using the MIT60 EoS (upper panel) and the 
LS EoS (lower panel). 
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FIG. 11: Direction and polarization averaged GW luminosity 
spectra for binaries with Mi = 1.2 Mq and M2 = 1.35 Mq 
in the case of the MIT60 EoS (black) and the LS EoS (red). 



becomes maximal. Note that wc shift the time axis to 
match i = at this point. 

The final plunge produces a rather complicated wave 
pattern. After it has taken place, one can see a quasiperi- 
odic signal from the ringdown of the postmerger rem- 
nant. This oscillation is reflected in the luminosity spec- 
tra. In Figs. [Til and [T31 we show the direction and 
polarization averaged luminosity spectra dEcw/df = 
2tt^ P {\h+\'^ + |/ixP) as computed from the Fourier 
transformed waveforms h for both polarizations. The 
very pronounced peak (or even multiple peaks) corre- 
sponds to the ringdown frequency of the postmerger rem- 



13 



^ -5 





/ mum 


MIT60 


v I 


/ ''IP' ""\ 

I K J ringdown 




1 — 

i 


1 1 1 

1 ^max 


— 1 

LS 




1/ ringdown 


/WWww 



-4 -2 2 4 6 8 

t [ms] 



FIG. 12: Gravitational-wave amplitudes for the plus polar- 
ization measured perpendicular to the orbital plane at a dis- 
tance of 20 Mpc for binary mergers with Mi = 1.35 Mq and 
M2 = 1.35 Mq, using the MIT60 EoS (upper panel) and the 
LS EoS (lower panel). 
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FIG. 13: Direction and polarization averaged GW luminosity 
spectra for binaries with Afi = 1.35 Mq and M2 = 1.35 A/q 
in the case of the MIT60 EoS (black) and the LS EoS (red). 



nant. This peak is located at a frequency /peak- Together 
with the maximum frequency during inspiral we now have 
at hand two characteristic features of the GW emission 
of compact object mergers. For a discussion and intro- 
duction of these quantities in the context of NS merg- 
ers, see where also the prospects of a detection and 
determination of these characteristics by GW detectors 
is briefly addressed. Here we will explore to what ex- 
tent these features can serve to distinguish NS mergers 
from SS mergers. It should be mentioned that the low- 
frequency part of our spectra (< 1 kHz) is not reliable 
because our simulations start only a few orbits before 



the plunge and therefore power in the low-frequency do- 
main from the preceding inspiral phase is missing. In the 
Newtonian limit the shape of the GW energy spectrum 
for the inspiral phase is given by a f~^^^ power law. 

Comparing the results for the binary with the 1.2 Mq 
and 1.35 Mq components one finds that the maximal 
frequencies /max that are reached during the inspiral are 
1.80 kHz for MIT60 and 1.57 kHz for LS (Fig.[Tni). The 
lower value of the NS model appears reasonable consid- 
ering the lower compactness of the initial stars compared 
to the quark stars. The ringdown signal after the merg- 
ing decays more slowly for the nuclear EoS (Fig.[Tni). We 
attribute this to the higher asymmetry of the merger rem- 
nant caused by the tidal disruption and suspect that the 
persistent oscillations may be supplied by the continued 
contraction of the remnant (Fig. [3]). Note also that the 
amplitude shortly after the merging is higher for the SS 
merger and one can clearly recognize a persistent mod- 
ulation of the signal. The origin of this low-frequency 
feature is unclear to us and should be examined by a 
mode analysis of the postmerger remnant. 

The spectra look qualitatively similar and in particular 
the peak frequencies of the postmerger ringdown are very 
close, both about 3 kHz (3.04 kHz for the NS binary and 
3.14 kHz for the SS binary). A remarkable difference is 
the pronounced gap below the peak frequency in the spec- 
trum of the case with the LS EoS. Also the reason for this 
feature should be investigated in detail by a mode anal- 
ysis. In general, a more detailed analysis of the spectra 
should reveal some more, less obvious differences. How- 
ever, this would require a detailed understanding of the 
mode excitation during the merging, which is beyond the 
scope of this work. It is also questionable whether such 
secondary differences can be of observational relevance 
for detectors available in the near future. 

In the case of symmetric configurations with two 
1.35 Mq stars the maximal frequencies during inspiral 
arc 1.92 kHz for MIT60 and 1.75 kHz for the nuclear LS 
EoS, and as before the GW amplitude shortly after the 
merging is higher for the MIT60 EoS, while the signal for 
the LS EoS is weaker but decreases more slowly (Fig. [T^ . 
A low-frequency modulation of the SS merger signal oc- 
curs similar to the waveform found for the asymmetric 
SS binary. As mentioned above, the merger remnant for 
MIT60 collapses after some revolutions and so the signal 
cannot be followed by our simulation any more. Nev- 
ertheless, the GW spectrum for the SS merger exhibits 
more power at high frequencies than the NS merger emis- 
sion and as in the asymmetric case the peak frequency 
of the SS merger is slightly higher. Again one recognizes 
a deep trough in the spectrum of the model with the LS 
EoS. Although present in all of our models with the LS 
EoS, it is unclear whether this is a universal feature and if 
it can be used to distinguish SS mergers from NS merg- 
ers in general (for the simulations employing the Shen 
EoS we also find a pronounced gap feature in the spectra 
except for two binary configurations; mergers described 
by the MIT80 EoS only form hypcrmassivc remnants for 
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FIG. 14: Maximal frequency during inspiral color-coded in 
kHz as a function of the total binary mass Mtot and of the 
mass ratio q for binary mergers with the MIT60 EoS. 
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FIG. 15: Maximal frequency during inspiral color-coded in 
kHz as a function of the total binary mass Mtot and of the 
mass ratio q for binary mergers with the MIT80 EoS. 



Mtot < 2.35 Mo). 

Since the binaries containing two stars with 1.8 Mq 
collapse promptly to BHs for both EoSs, one can only 
consider the maximal inspiral frequency as character- 
istic quantity. For the LS EoS one finds /max ~ 
2.16 kHz, while for the SS binary the maximal frequency 
is 2.20 kHz. Since the stars of this binary setup have the 
same radius for both EoSs, the similarity of these values is 
not unexpected. In [2l| it was suggested that the mass- 
radius relation and thus the EoS of compact stars can 
be determined from deviations of the GW signal from 
the point-particle behavior. However, this proposal re- 
lies on the close relation between the signal properties 
and just the stellar radius as the crucial stellar param- 
eter, and neglects the potential influence of a different 
inner structure of compact stars with the same mass and 
radius. The same values of /max, one of the most remark- 
able features of the inspiral phase of our I.8M0-I-I.8M0 
models, confirms that the inner structure has indeed a 
minor influence even on the dynamics of the late stages 
of the inspiral phase. Also in our study we find that 
the mass-radius relation is the EoS property that most 
sensitively determines /max, because for an NS merger 
with the stiffer Shen EoS, where a 1.8 Mq NS has a 5 
km bigger radius, we obtain /max = 1-62 kHz. On the 
other hand, a disadvantage of the insensitivity to the 
inner structure may be that it requires several "radius 
measurements" to constrain the EoS, because one detec- 
tion may not provide unambiguous information about the 
EoS. 



B. Binary parameter dependence 

Figures [T31 HH HZ] and (TS] present in overview the 
characteristic frequencies for all models that we com- 
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FIG. 16: Maximal frequencies during the inspiral phase as a 
function of the compactness of the initial stars for all config- 
urations with a mass ratio of g = 1 and all EoSs used in this 
study. The compactness is defined as the gravitational mass of 
a single star in isolation divided by its radius in Schwarzschild 
coordinates. 



puted with the MIT60 and MIT80 EoSs. The results 
are displayed in dependence of the total binary mass 
Mtot = Ml -\- AI2 and of the mass ratio q = Mi/M2- 
Note that the peak frequency of the postmerger ring- 
down is given only for the configurations that form a 
hypermassive remnant. 

For the maximum frequency of the inspiral we find a 
monotonic, but not exclusive, dependence on the total 
mass of the binary as shown in Figs. [T3] and [151 This 
trend is also present when considering /max as a function 
of the compactness of the initial stars as given for all 
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FIG. 17: Peak frequency of the postmerger ringdown color- 
coded in kHz as a function of the total binary mass Mtot and 
of the mass ratio q for binary mergers with the MIT60 EoS. 
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coded in kHz as a function of the total binary mass Mtot and 
of the mass ratio q for binary mergers with the MIT80 EoS. 



symmetric configurations in Fig. 1161 The more compact 
the stars are, the higher the frequency is that is reached 
during the inspiral. Remarkably, /max seems to follow an 
approximately linear function independent of the EoS. 
This confirms that the inner structure has a small effect 
on the inspiral phase. Only the relatively massive con- 
figurations describe a kink that deviates from the nearly 
linear behavior. The mass ratio also seems to have an in- 
fiuence, because in Figs.lHlandfTSl f^,^-, slightly decreases 
for lower values of q. Qualitativel y th e same result was 
found for the mergers of NSs (see [23l ) . Our findings for 
the maximum GW frequency for the inspiral phase of SSs 
are also in agreement with [25| , where symmetric binaries 
were considered and a similar behavior was seen. 



FIG. 19: Peak frequencies of the postmerger ringdown as a 
function of the remnant compactness for all EoSs considered 
in this study. While the surface of SS merger remnants can be 
clearly determined, we define the size of NS merger remnants 
arbitrarily by the isodensity surface with p = 10^^ g/cm"^. 
The luass of the remnant is approximated by the total mass 
Mtot ~ Ml + M2- The compactness is then given as A/tot 
divided by the remnant radius in the equatorial plane in 
isotropic coordinates at the end of the simulation when an 
axisymmetric hypermassive object has formed. Note that the 
differentially rotating remnant is highly deformed with an axis 
ratio of 2:2:1. For each EoS an approximately linear relation 
between /peak and the compactness is obtained. The outliers 
are relatively massive configurations, which collapse to a BH 
during the simulation time. 



Figures [T7] and [TS] reveal that for both EoSs the fre- 
quency of the postmerger peak increases monotonically 
with the total mass of the system. We interpret this de- 
pendence as a consequence of the growing compactness 
of the hypermassive remnant (see Fig. [T5)) , which is de- 
termined by the total system mass. The compactness in 
turn determines the rotational period and the frequency 
of oscillations of the deformed object that produces the 
GW signal. Corresponding results for neutron mergers 
have been reported by the authors of Ref. [l^], whose 
findings we confirm by our additional NS merger simula- 
tions. 

In addition, one observes that the mass ratio hardly 
affects the value of /peak- As mentioned above, the dy- 
namics of the postmerger remnant are determined by the 
total mass of the single object that forms in the merger 
and is therefore nearly independent of the binary mass 
ratio. 



Equation of state dependence and thermal 
effects 



Obviously, a comparison of the influence of the differ- 
ent EoSs for SQM and for ordinary nuclear matter is of 
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great importance. Figure [20] compiles the characteristic 
frequencies for the MIT60, MIT80, LS, and Shen models. 
The SS EoSs exhibit the tendency of being located in the 
upper right part of the plot, while LS and Shen give rel- 
atively low /max between about 1.1 kHz and roughly 1.8 
kHz for the maximum values. One should keep in mind 
that wc did not compute the same sample of binary con- 
figurations for all EoSs, because the different EoSs sup- 
port differently high masses and determine the system 
properties for which a hypermassive remnant forms and 
/peak can be extracted. However, a comparison of the 
same binary setup for all EoSs is possible in the cases of 
the colored symbols in Fig. [20] (blue: 1.35Af0-hl.35Mo 
configuration; red: 1.2M0-I-1.2M0 configuration) and 
the trends visible arc compatible with the rest of the 
data. 

Again the diagram can be understood by means of the 
compactness of the merging stars and postmergcr rem- 
nants: The more compact the stars or the merger rem- 
nants are, the higher are /max and /peak, respectively (see 
Figs. [TBI and [T^ and compare also with Fig. [5]). Remark- 
ably, MIT60 and LS show an overlap of /max and /peak 
values, which means that they can not be distinguished 
by a determination of these characteristic frequencies in 
GW measurements. 

In this context we refer to the discussion in (56| , which 
predicts an enhanced observable flux of strangelets in cos- 
mic rays if the strange matter hypothesis was true and 
quark matter was described by an EoS similar to MIT60. 
In contrast, if the properties of SQM caused SSs to be 
as compact as those obtained with the MIT80 EoS, this 
should be clearly visible in the high characteristic fre- 
quencies of the GW signal. On the other hand, a nuclear 
EoS like the Shen EoS appears to be distinguishable from 
SQM because the MIT60 EoS represents a lower limit 
with respect to the compactness of SSs. 

The existence of strangelets may be also crucial to dis- 
cern self-bound SSs from hybrid stars, i.e. from NSs with 
cores of quark matter that is not absolutely stable and 
thus cannot form self-bound quark matter nuggets (se e 
e.g. 0). The rich diversity of hybrid stars (see [66j ) 
also includes objects with stellar parameters comparable 
to those obtained with the LS EoS or the MIT60 EoS. 
Therefore, hybrid stars might populate the "LS-MIT60- 
area" in Fig. 1201 However, for many hybrid star models 
we expect also the GW signals to be discriminating, since 
these objects can have a compactness, as a characteristic 
quantity determining the GW features, well below the 
one of NSs described by the LS EoS [6l|. This clearly 
distinguishes them from SSs. For an investigation of the 
differences of the GW signals of ordinary NSs and hybrid 
stars we refer the reader to [s^l . 

The comparison of the GW signals from cold and warm 
simulations confirms the findings of the previous section 
that the additional pressure support by nonzero temper- 
ature effects is important for NS mergers, while for the 
SS models the change of the gravitational mass in the 
zero temperature simulations is the most relevant aspect 
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FIG. 20: Maximal frequencies during inspiral and peak fre- 
quencies of the ringdown of the postmerger remnant for all 
configurations that form a hypermassive object. Results for 
all EoSs used in this study are plotted. The blue symbols 
correspond to the binary configurations with two 1.35 Mq 
stars, the red ones mark the frequencies for merger events of 
two stars with 1.2 Mq. Note that there is no blue symbol for 
the MIT80 EoS because the configuration with two 1.35 Mq 
stars leads to an immediate BH formation. The compilation 
includes data from [221. 
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FIG. 21: Enclosed gravitational mass within an ellipsoid 
with the semiaxes a — x, b = x and c — x/2 for the 
1.2Af0-|-1.35MQ binaries. Shown are results with the LS 
EoS and the MIT60 EoS including (red curves) and neglect- 
ing (black curves) thermal efi^ects. Note that within general 
relativity the gravitational mass is only defined in isolation. 
Here Mgrav denotes the enclosed contribution to the Arnowitt- 
Deser-Misner mass, neglecting the extrinsic curvature terms. 
Using the enclosed rest mass yields a similar result. The dis- 
tance x is given in isotropic coordinates. 
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(see Tabic H]). As one might expect, /max is very simi- 
lar because thermal effects do not play a role during the 
inspiral. In fact, the variation remains below the rela- 
tively large uncertainties of about 0.1 kHz for /max that 
are connected with the numerical determination of the 
merging time. 

Comparing the ring-down signals, one finds that for 
the models with the LS EoS /peak is higher in the cold 
than in the warm simulations (Table |I| . Figure [5T] pro- 
vides the explanation of this result. It displays the mass 
enclosed by ellipsoidal surfaces with varied distance from 
the center of the remnant. One can see that for the nu- 
clear LS EoS the cold remnant is more compact because 
of the missing thermal pressure support. Consequently, 
the peak frequencies are higher. 

In contrast to the NS models, we find lower peak fre- 
quencies in the cold SS merger simulations. As visible 
in Fig. [21] in the case of SSs the cold remnants are less 
compact due to the lower gravitational mass, while the 
remnant size is essentially unaltered since the thermal 
pressure support is not important for SQM. The reduced 
compactness explains the lower peak frequencies for SS 
merger simulations with T = 0. 



D. Occurrence of a prompt collapse 

The occurrence of a prompt collapse might be another 
feature that can be used to discern SSs from NSs. From 
GW measurements the chirp mass and most probably 
the individual masses of the binary components can be 
deduced, and therefore the total mass and the mass ra- 
tio. In addition, by analyzing the power at higher fre- 
quencies it should be possible to judge from a measure- 
ment whether a hypermassive postmerger remnant or a 
BH formed. This idea has already been brought forward 
in [g^l. Therefore we checked how the mass limit for a 
prompt collapse depends on the EoS. For the sake of clar- 
ity we consider only equal-mass binaries, but we expect 
the analysis to also hold for asymmetric systems. 

For the LS EoS the minimum total mass A/coiiapsc for 
configurations that undergo a prompt collapse instead of 
forming a hypermassive object is between 2.8 Mq and 
2.9 Mq. The corresponding mass limit for the MIT60 
EoS is slightly higher than 2.7 Mq, but below 2.8 Mq. 
As mentioned above, the symmetric binary with Mtot = 
2.7 Mq does not collapse promptly but only ~ 4 ms 
after the hypermassive object has formed. Therefore, 
by this property a marginal discrimination of these two 
EoSs is possible. For MIT80 we find McoUapsc < 2.4 Mq, 
which distinguishes this EoS clearly from the majority 
of nuclear EoSs, because the LS EoS represents a fairly 
extreme case in terms of softness, NS compactness and 
McoUapsc- The threshold mass for the Shen EoS is about 
3.5 Mq. For asymmetric systems we refer to Fig. 1 in 
[56j . which shows the borderline between binary config- 
urations that result in a prompt collapse and those that 
lead to the formation of a hypermassive object. 



The criterion discussed here can be considered as very 
simple and straightforward in its application to observa- 
tional data. This makes it attractive for an analysis of 
merger signals and for the determination of basic con- 
straints on the EoS. A drawback, however, may be the 
fact that GW measurements for a larger set of events will 
be needed to constrain the mass limit for the direct BH 
formation. Alternatively, the delay time between plunge, 
defined as the time when the GW amplitude becomes 
maximal, and collapse, which ultimately occurs for ev- 
ery hypermassive object, will be discussed as source of 
information in [ssf . For the models described in detail in 
Sect. HVl we report the delay times Tdoiay in Table |T] for 
the cases where the collapse occurs during the simula- 
tion. Otherwise the time until the end of the simulation 
gives an lower limit on Tdoiay, but the true value may be 
significantly larger. 



E. Gravitational-wave luminosity 

As argued above it is not possible in all cases to dis- 
tinguish NS mergers from SS mergers by extracting rela- 
tively simple characteristic frequencies from a GW mea- 
surement. If the mass-radius relations of compact stars 
were similar even only in the most relevant range of stel- 
lar masses, the frequencies /max and /peak niay be in the 
same range of values for different EoSs. A more detailed 
analysis may yield additional discriminating information 
but will also require a higher quality of the signal deter- 
mination. 

Besides the method of matched filtering, which re- 
quires the knowledge of the signal searched for, analy- 
sis pipelines of GW detectors have the ability to mea- 
sure excess power in time- frequency domains josl [69| . 
which allows to search for unmodeled signals. Ratios 
of frequency-integrated energies in the luminosity spec- 
trum, for example, can be considered as signal character- 
istics that can be extracted from the observational data 
without knowing the signals in detail. For instance one 
could consider the ratio between the energy Aii'pm ra- 
diated away by the postmerger remnant within the first 
5 ms after merging and the energy A_Ein emitted in the 
last 3 ms of the inspiral before the coalescence. The time 
of merging is defined as the time when the GW ampli- 
tude becomes maximal. This quantity was introduced 
for characterizing binary NS merger GW signals in [2^ . 
The energies can be determined either by an integral over 
the luminosity spectrum or more easily by a time inte- 
gral over the luminosity given directly by the quadrupole 
evolution dE/dt = l/5(Qjj(5jj). Both ways of evaluation 
yield quantitative agreement within 20% when applied 
to the results of our numerical simulations. For simplicity 
we consider the latter expression. 

Figure [22] shows the ratio of AiJpm to AiJin as a func- 
tion of the peak frequency /peak of the postmerger ring- 
down for all possible binary configurations of stars with 
1.2 Mq and 1.35 Mq, and for the EoSs of Shen, LS and 
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FIG. 22; Ratios of the GW energy emitted by the postmerger 
remnant within 5 ms after the coalescence to the energy ra- 
diated during the last 3 ms of the inspiral as a function of 
the peak GW frequency of the postmerger remnant. Shown 
are results for the binary configurations with 1.2 Mq and 
1.2 Mq (left datapoints of each EoS cluster), with 1.2 Mq and 
1.35 Mq (central datapoints), and with 1.35 Mq and 1.35 Mq 
(right datapoints) for the Shen, LS, and MIT60 EoS. 



MIT60. Note that /peak increases with the total mass of 
the system as shown in Figs. [171 and[T51 allowing for an 
unambiguous identification of the three configurations in 
this plot (the least massive binaries are located to the 
left of the EoS-clusters of datapoints). Clearly one can 
see that the ratio Ai?p,„/AiJin is a valuable measure to 
discriminate NS mergers described by the LS EoS from 
SS mergers with the MIT60 EoS. One may suspect such 
differences already from a comparison of the wave ampli- 
tudes during the postmcrging phase plotted in Figs. 1101 
and[ll 

Motivated by the differences in the time evolution of 
the GW amplitudes visible in Figs. [TU] and [T^l we also 
computed the energy ratios for postmerging time inter- 
vals Atpm of 3 ms, 5 ms, and 10 ms. Figure [231 reveals 
a noticeable discrepancy between SS mergers and coales- 
cence events of ordinary NSs. Although the energy ratio 
is initially higher for SS postmerging remnants, the lumi- 
nosity of these objects decays much faster and the value 
of AEpn-i/ AEin saturates after about 5 ms. In contrast, 
this ratio increases for NS mergers over emission times 
even longer than 5 ms after the coalescence. Therefore, 
also the decay time of the postmerging emission should 
be considered as a feature that can be used to discern the 
GW signal of SS mergers and NS mergers, in particular in 
cases when the characteristic signal frequencies /max and 
/peak do not allow for an unambiguous discrimination. 

As visible in Figs. [TTl and [T3l the GW luminosity spec- 
tra of NS mergers employing the LS EoS show a promi- 
nent gap at frequencies slightly lower than /peak. This 
feature occurs in all of our models with the LS EoS, while 



FIG. 23: Ratios AEpm/ AEin for different emission periods 
Atpm of GWs during the postmerging phase. 



it is less pronounced in the cases of mergers with the 
MIT60 EoS. The gap is also visible in the Fourier trans- 
formed waveform and it occurs for each emission direc- 
tion. This suggests the use of the gap to break the de- 
generacy of the LS and MIT60 EoSs in /^ax and /peak 
(Fig. [20|) ; we focus on these two EoSs because the char- 
acteristic frequencies do not allow their discrimination as 
in the cases of the other EoSs. 

We introduce frequency intervals with a given width 
A/ placed centered around /peak and /gap. The latter 
frequency is defined by the position of the local mini- 
mum that defines the gap in the luminosity spectrum 
(Figs. [TTl and [T5)) . The values of /gap for the more mas- 
sive binary configurations are listed in Table HI while for 
the 1.2Mq+1.2Mq binary /gap equals 2.37 kHz for the 
LS EoS and 2.57 kHz in the case of the MIT60 EoS. In- 
tegrating the luminosity spectra in these frequency win- 
dows and comparing the energies provides a characteris- 
tic quantity for a quantitative description of the gap fea- 
ture of NS merger models with the LS EoS. This energy 
ratio Ai^peak/AE'gap may be deduced easily from obser- 
vational data by analyzing the excess power in different 
frequency domains [68l . |69| . Figure [Ml shows the ratio 
Ai?pcak/Ai?gap of binary mergers with the LS and MIT60 
EoSs for different choices of A/. Equal colors correspond 
to the same value of A/. For smaller frequency windows 
the difference in the ratios becomes more enhanced as 
one expects from Figs. [TTl and [121 One can clearly see 
that for all choices of A/ the ratio AE'pcak/AE'gap allows 
to distinguish the GW emission of mergers described by 
the LS EoS from the signal of the SS coalescence events 
with the MIT60 EoS. 
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FIG. 24: Ratios ASpeak/ASgap of the GW energy emitted 
in a frequency window A/ around the peak frequency to the 
energy emitted in a frequency window A/ around the gap 
frequency as a function of the peak frequencies of the post- 
merger ringdown. The colors correspond to differently chosen 
frequency windows A/ (see main text), therefore lines of the 
same color are to be compared (blue: 250 Hz; red: 100 Hz; 
black: 50 Hz). Note the logarithmic scale of the ratio. Shown 
are results for the binary configurations with 1.2 Mq and 
1.2 Mq (left datapoints), with 1.2 Mq and 1.35 Mq (cen- 
tral datapoints), and with 1.35 AIq and 1.35 Mq (right dat- 
apoints) for the LS EoS (dotted hues) and MIT60 EoS (sohd 
lines) . 



VI. SUMMARY AND CONCLUSIONS 

We have reported on simulations of SS mergers and 
discussed our results in comparison to NS coalescence 
events. In particular, we focussed on the GW signals and 
their potential to provide information about the existence 
of SSs. 

We found that the dynamical behavior of SS mergers is 
fundamentally different, which can be understood by the 
higher compactness of SSs and their merger remnants, 
the self-binding of SQM and the different influence of 
thermal effects. While NS mergers form a dilute halo or 
toruslike structure around a dense, hypermassive, differ- 
entially rotating remnant, the remnant of merging SSs 
is bounded by a sharp surface as were the initial stars. 
Only by the formation of thin tidal arms relatively late 
in the evolution of the remnant, matter gets shed off the 
central object and forms a fragmented thin disk in the 
equatorial plane. 

In order to estimate the importance of thermal ef- 
fects during NS and SS mergers we compared our models 
with simulations where zero temperature was imposed 
by extracting the thermal energy as most efficient cool- 
ing would do. For both kinds of stars we observed a 
sensitive dependence of the dynamics of the system on 
thermal effects, affecting for instance the development 



and structure of the outer remnant parts of NS merg- 
ers. While the estimated relic torus mass after collapse 
of the remnant to a BH is in general higher if thermal ef- 
fects were neglected, a nonzero temperature of NS matter 
leads to an inflated, dilute halolike torus in contrast to 
a much thinner, more disklike structure for T = 0. For 
SS mergers we found a shorter time delay for the BH 
formation, whereas in the case of NSs this time interval 
is stretched as reported in [65|. This difference can be 
understood as a consequence of the additional gravitat- 
ing effect of the thermal energy, which in the case of very 
dense SSs is not overcompensated by thermal pressure 
effects (different from the NS case). 

The analysis of the GW signals emitted by SS merg- 
ers revealed that already by means of relatively simple 
characteristic features of the signal it may be possible to 
decide whether SS or NS mergers produced the emission. 
In particular, we found that the maximal frequency dur- 
ing the inspiral and the frequency of the ringdown of the 
postmerging remnant arc in general higher in the case 
of SS mergers in comparison to NS mergers. Whether 
this criterion can be used to finally decide on the strange 
matter hypothesis depends on the particular stellar prop- 
erties associated with the EoS of high-density matter. 
For a similar mass-radius relation within a certain mass 
range, meaning relatively compact NSs or less compact 
SSs (the LS-MIT60 scenario), the determination of these 
frequencies might not be decisive. However, taking ad- 
ditional characteristics of the GW luminosity into con- 
sideration will allow for a discrimination. In this con- 
text we discussed the occurrence of a prompt collapse, 
the ratios of the GW energy emitted in the postmerger 
phase to the energy radiated away during the inspiral 
phase, the growth rate of the energy emission associated 
with the postmerger ringdown signal, and the appear- 
ance of a gap in the GW luminosity spectrum. In ad- 
dition, it was shown in [56j that cosmic ray experiments 
will yield information about the strange matter hypoth- 
esis and may clarify the situation, especially in the case 
when GW signals are not conclusive as in the LS-MIT60 
scenario. Therefore we expect that despite the rather un- 
certain event rate the upcoming advanced GW detectors 
LIGO and VIRGO may provide valuable data to decide 
about this long-standing question, and it appears likely 
that one will thus gain fundamental information about 
the properties of high-density matter. Moreover, future 
GW detectors like the Einstein telescope (T^I and the 
DUAL detector [7l| will have a higher sensitivity in the 
high-frequency domain above 1 kHz, where characteris- 
tic features occur. Signals measured with these future 
instruments will therefore allow for a more detailed anal- 
ysis. 

Finally, we would like to point out that also other forms 
of self-bound matter have been discussed in the literature 
(see e.g. [U, [zl, [l^ ) ■ For instance pion condensates may 
lead to stellar objects similar to SSs, and also abnor- 
mal nuclei similar to strangelets may exist in this case. 
Therefore, the study presented in this work should be 
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generalized to these states of matter, where quantitative 
considerations require the knowledge of the specific EoS. 
However, for stellar properties comparable to those of 
SSs we expect qualitatively the same behavior and the 
same GW features as for SSs. From this discussion it 
is clear that only a multimessenger approach can bring 
a final answer. In addition to GW measurements and 
cosmic ray experiments, this may in particular include 
the observation of isolated compact stars to derive the 
cooling history of these objects, which is significantly dif- 
ferent for NSs, SSs and self-bound pion-condensed stars 
(see [53,[2HZ3)- 

Several improvements and extensions may supplement 
this first study of SS mergers in the future. Besides re- 
finements in the methodical approach like a fully rela- 
tivistic treatment and the inclusion of magnetic fields 
(only important during the merging phase if very strong, 
but must be expected to grow afterwards and affect the 
ringdown phase) and radiative effects, one would also 
like to account for the effects of quark interaction and 
color superconductivity, employ models beyond the MIT 
bag model, and discuss other forms of self-bound mat- 
ter. Moreover, a determination of the modes excited in 



the merger remnants may be a promising way to develop 
a better understanding of the origin of the characteris- 
tic properties of the GW spectra. Furthermore a study 
of SS-BH mergers will yield further insights about the 
behavior of such systems in contrast to NS-BH binaries. 
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